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The top-quark mass, nit, is an intrinsic parameter of 
the standard model (SM) of particle physics and is of par- 
ticular importance due to its strikingly large value. As a 
result, the top quark has a large effect on radiative correc- 
tions to electroweak processes and has a Yukawa coupling 
to the Higgs field of 0(1), which may provide insight into 
the mechanism of electroweak symmetry breaking [T]. 

The Higgs boson mass, win, is not predicted by the 
SM, but constraints on its value can be derived from 
the calculation of radiative corrections to the W boson 
mass, niw, and from the values of other precision elec- 
troweak variables [2J. These corrections depend primar- 
ily on lnm# and m 2 , and thus precision measurements 
of mw and m t provide important constraints on nin . 

The dominant top-quark production process is pair 
production via the strong interaction. At Fermilab's 
Tevatron this process is initiated by pp collisions at 
center-of-mass energy yfs — 1.96 TeV. Because of its 
large mass, the top quark decays rapidly with lifetime 
r t ~ 10 -25 s [3] — fast enough that it has essentially no 
time to interact and may be considered as a free quark. 
This allows a direct measurement of its mass from the 
daughter particles from its decay, and as a result nit 
has the lowest relative uncertainty of all of the quark 
masses [I]. 

In the SM top quarks decay via the weak inter- 
action, predominantly to W bosons and b quarks as 
ti — > W + b W~b. W bosons decay into lower-mass 
fermion-antifermion pairs: a charged lepton and a neu- 
trino (W + — > Ivt or W~ — > Ivg), "leptonic decay"; or 
an up-type quark and a down-type quark (W + —> qq 1 
or W~ — > qq'), "hadronic decay". The result presented 
here uses the lepton+jets decay channel (with qq'bii'ib or 
Ivibqq'b in the final state), where one of the two W bosons 
decays leptonically into an electron or a muon, and the 
other decays hadronically. All the quarks in the final 
state evolve into jets of hadrons. Events with tau lep- 
tons are not selected directly, but may contribute a few 
percent of the total sample via leptonic cascade decays 
or fake jets. The most recent m t measurements obtained 
at the Tevatron using the lepton+jets topology are re- 
ported in Ref . !5| , while the results of an earlier version 
of the present analysis using 955 pb -1 of integrated lu- 
minosity are reported in Ref. [BJ. The distinctive feature 
of this analysis is the use of matrix element calculations 
to describe the dominant background contribution. The 
result presented here uses a more than three times larger 
data sample than the earlier version, and employs a more 
detailed likelihood function. 

The leptons and jets resulting from the top-antitop 
quark pair (ti) decay are detected in the CDF II general- 
purpose particle detector that is described in detail else- 
where [7] . Azimuthally and forward-backward symmetric 
about the beam-line, the detector contains a high preci- 
sion particle tracking system immersed in a 1.4 T mag- 
netic field and surrounded by calorimetry, with muon 
detectors on the outside. A right-handed spherical co- 
ordinate system is employed, with the polar angle 9 mea- 



Electron E T > 20 GeV |t?| < 1.1 
or Muon p T > 20 GeV/c jjjj < 1.0 
E T Et > 20 GeV jjjj < 3.6 

Jets E T > 20 GeV \r)\ < 2.0 
Four jets; at least one from a b quark 



TABLE I: Event selection criteria. 

sured from the proton beam direction, the azimuthal an- 
gle (f> in the plane perpendicular to the beam-line, and 
the distance r from the center of the detector. Trans- 
verse energy and momentum are defined as Et = E sin 9 
and pt = psin0, where E and p denote energy and mo- 
mentum. Pseudorapidity is defined as n = — In tan (9/2). 

This measurement makes use of CDF II data col- 
lected between February 2002 and August 2008, repre- 
senting approximately 3.2 fb _1 of integrated luminos- 
ity. The event selection criteria (Table [I]) are tuned 
to select the lepton+jets final-state particles, requiring 
that each event must have exactly one high-i^y elec- 
tron or high-px muon, exactly four high-I^r jets, and 
a significant amount of missing Et, /Et [1], charac- 
teristic of the undetected neutrino. Jets are recon- 
structed using a cone algorithm [5] , with the cone radius 

AR = \J{Anf + (A0) 2 = 0.4. At least one of the four 
jets must be identified as originating from a b quark via 
the SECVTX algorithm [TU], which detects displaced 
secondary vertices characteristic of the decay of long- 
lived b hadrons. A total of 578 events are selected, of 
which 76% are expected to be ti events (Table [TTJ) . Of 
the 24% of events expected to be background, it is pre- 
dicted that 69% arise from the production of a W bo- 
son in conjunction with 4 jets (W+jets), 19% come from 
multi-jet QCD production (non-14 7 ), while the remain- 
ing 12% are from sources such as diboson and single- 
top-quark production. These fractions are estimated us- 
ing theoretical cross-sections, Monte Carlo (MC) sim- 
ulated events, and data. The ti events are generated 
using the Lund Monte Carlo program PYTHIA [llj . 
with a top-quark mass of 175 GeV/c 2 and a ti produc- 
tion cross section of 6.7 ± 0.8 pb [12]. The W+jets and 
Z+jets events are generated using the ALP GEN gen- 
erator jT3] while the single-top-quark events are gener- 
ated using the MADEVENT package [14], in both cases 
also using PYTHIA to perform the parton showering 
and hadronization. Diboson events are also generated 
using PYTHIA. In addition, data are used for non-11 7 
events [15] . 

This analysis employs an unbinned maximum likeli- 
hood method [THl E?J . The mt-dependent probability 
density function (p.d.f.) is calculated for each event in 
the data sample: 

V(k) = v sig P s (k) + (1 - v sig )P b (k), (1) 

where k = (Ei,pi) represents the measured kinematic 
quantities of the event, P s and Pb are respectively the 
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sample # of events % of total % of bkg 



tt signal 


425.0 ± 58.9 


76.0% 


- 


VK+jets 


92.6 ± 15.9 


16.6% 


69.0% 


non-W 7 


25.0 ± 12.5 


4.5% 


18.7% 


single top quark 6.6 ± 0.4 


1.2% 


4.9% 


diboson 


6.0 ± 0.6 


1.1% 


4.5% 


Z+jets 


3.9 ± 0.5 


0.7% 


2.9% 


total 


559.2 ± 67.0 


100% 





Observed 578 



TABLE II: Number of expected signal and background events, 
corresponding to the total integrated luminosity of 3.2 fb _1 . 
The percentages are used when generating Monte Carlo sim- 
ulated experiments. 

normalized p.d.f.s for signal and background events, and 
f s ig is the signal fraction parameter (constrained < 
^sig < 1)- Signal events are defined as events consis- 
tent with qq — > tt production and tt decay into the lep- 
ton+jets channel, as described by the leading-order (LO) 
matrix element evaluated by Mahlon and Parke [IE] . 
Background events are assumed to be described by a ma- 
trix element for ly+jets production, which is calculated 
using a sum of 1286 VF+4-partons amplitudes for 592 
subprocesses encoded in the VECBOS MC event gener- 
ator [19] . This approximation does mean that there are 
some events that, in principle, are not described by either 
P s or Pt, including non-W, single top, diboson, Z+jets, 
and W + &6+2-partons events, as well as W+jets events 
from VF+0-, 1-, 2- and 3-partons processes. However, 
studies with MC simulated events show that the ratio 
Pb/P s calculated for all of these event types is similar 
to that for VF+4-partons events, and that, in practice, 
such events mostly contribute to the likelihood function 
via the P^ term and do not add any more bias than the 
W^+4-partons events or than the poorly reconstructed tt 
events themselves [20] . Any residual bias in the measured 
top-quark mass is removed at the end, as described later 
in the paper. 

The signal and background p.d.f.s, P s and Pi,, are 
constructed in analogous fashions, starting with the 
appropriately normalized parton-level differential cross- 
section [I], da s or dab, which is then convolved with 
parton distribution functions (PDFs) and a jet-to-parton 
transfer function W(k,>c). P s is thus given by 

1 njP 1 1 

P s (k;m t , Aj E s) = V —, r-r-, T rX 

(2) 

J da s {x;m t ) dx^dx^ W(k, k; A JE s) f(x^)(x^), 

where n = (ej,7Tj) represents the actual event parton- 
level kinematic quantities corresponding to the measured 
quantities k, and parameter Ajes is defined in the next 
paragraph. The PDFs /(^Bj) define the probability den- 



sity for a colliding parton to carry a longitudinal momen- 
tum fraction XBj and are given by CTEQ5L [21] , A s is 
the mean acceptance function for signal events, a normal- 
ization term that is the consequence of the constriction 
of the phase-space of the integral by the event selection 
cuts and by the detector acceptance. The average over 
the jet permutations, rij p , is due to ambiguity in assign- 
ing final state jets to partons. The fact that the two light 
quarks in the final state are indistinguishable allows the 
reduction from the original 24 permutations to 12 in the 
expression for P s , and the 6-tagging information allows 
a further reduction to 6 assignments for events with one 
identified 6-jet and 2 for events with both 6-jets identi- 
fied. In the similar expression for P^, all 24 permutations 
are averaged. 

The jet-to-parton transfer function W{k, x) is a p.d.f. 
describing the probability density for an event with out- 
going partons and charged lepton with x to be measured 
as reconstructed k. The charged lepton is assumed to 
be well-measured, allowing the use of a Dirac ^-function 
to represent the mapping between its parton-level mo- 
mentum, 7?£, and its reconstructed momentum, p£. For 
the four jets, the function is obtained by parameteriz- 
ing the jet-to-parton mapping observed in fully simulated 
PYTHIA tt events. These events contain all of the in- 
formation about the original partons as well as the mea- 
sured jets. The simulation includes physical effects, such 
as radiation and hadronization, as well as the effects of 
measurement resolution and of the jet reconstruction al- 
gorithm. The parameterization is made in two parts that 
are assumed to be independent: the energy transfer func- 
tion We, describing the jet energies E, and the angular 
transfer function Wa, describing the mapping for the jet 
angles. The jet-to-parton transfer function is thus given 
by 

W(k,x;A JES )^5 3 (p e ~7f e )W A (3) 

The reconstructed jet energies, Ei, used in the func- 
tion We are not just the raw calorimeter energy deposits, 
but are first calibrated so that they represent the com- 
bined energies released in the calorimeter by the many 
particles constituting each jet. This is achieved using 
the CDF jet energy scale (JES) calibration (55], which is 
subject to a significant systematic uncertainty. The un- 
certainties of individual jet energy measurements, a(Ei), 
are therefore correlated, and their fractional JES uncer- 
tainty, a(Ei)/Ei, is typically ~ 3%. If this were included 
as a systematic uncertainty on the measured m t it would 
reduce the measurement precision drastically; in fact, 
each 1% of fractional JES uncertainty would add about 
1 GeV/c 2 uncertainty to the measured m t [23]. How- 
ever, such a treatment overestimates the uncertainty be- 
cause the energies of the two daughter jets of the hadroni- 
cally decaying W boson can be constrained based on the 
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FIG. 1: The reconstructed 2-jet invariant mass of the hadron- 
ically decaying W boson, mw, for measured jet angles (solid 
line) and for parton- level angles (dotted line), obtained after 
assuming the primary parton energy as jet energy. For ease 
of comparison, the parton-level distribution is normalized so 
that the maxima of the two distributions are the same. 



known W boson mass. Applying this constraint to all 
events in the data sample while allowing the jet energies 
to be shifted results in the in situ measurement of the 
JES correction, A JES , denned as the number of a(Ei) 
values by which the energy of each jet is shifted in the 
likelihood fit. This effectively re-calibrates the measured 
jet energies based on the known W boson mass and re- 
places a large component of the JES systematic uncer- 
tainty with a much smaller statistical uncertainty on the 
A jes. The A jes dependence of the jet energies is in- 
cluded in the parameterization of the function We- This 
parameterization is made in eight bins in pseudorapid- 
ity \r)\, separately for light and 6-jets, using a sum of two 
Gaussians as a function of the difference between the par- 
ton energies and the corrected jet energies as measured 
in a sample of PYTHIA it events that pass the same 
selection criteria as the data. 

In an earlier version of this analysis [5J, the jet-to- 
parton transfer functions for all jet angles were approx- 
imated by Dirac (^-functions. The introduction of the 
function Wa was motivated by a discrepancy noticed in 
simulated tt events in the 2-jet effective invariant mass 
of the hadronically decaying W boson, mw- Even when 
the true simulated parton-level jet energies are used, in- 
stead of the corresponding reconstructed detector-level 
values, the use of the measured jet angles rather than 
their parton-level values causes a significant shift of the 
reconstructed mw from its nominal value, as illustrated 
in Fig. [TJ 

There is also a negative skewness in the distribution 
for measured angles, and since parton-level jet energies 



are used, the observed effects are due to the differences 
between the measured angles and the parton-level angles 
alone. The peak of the mw distribution, when fit by 
a Breit-Wigner distribution, corresponds to a W boson 
pole mass of 79.5 GeV/c 2 , a -0.9 GeV/c 2 shift from its 
parton-level value of 80.4 GeV/c 2 . This is found to be 
a result of a correlation between the measured jet direc- 
tions: the measured angle, a\2, between the two jets is, 
on average, reduced so that the two jets appear closer 
together than their parent partons, which can be seen 
in Fig. [2] Since the apparent W boson mass is utilized 
to measure Ajes and thus calibrate the measured jet 
energies, a jet-to-parton transfer function describing the 
change in the angle a.\2 is important in making an accu- 
rate measurement of Ajes an d thus the top-quark mass. 
The function Wa also describes a much smaller correla- 
tion effect seen in the angle awb between the hadronic- 
side 6-jet and the hadronically decaying W boson. The 
function Wa is thus parameterized using two different 
functions, W\ 2 and W A vb ', describing the mappings for 
the angles a,\i and awb- The remaining angles describe 
resolution effects rather than the correlations and, due to 
computational constraints, are assumed to be well mea- 
sured with their contributions to Wa approximated by 
Dirac delta functions. 

The functions W^ 2 and W\ vh are both fit using a sum 
of a skew-Cauchy distribution and two Gaussians, de- 
scribing the change in the cosine of the relevant angle, 
A cos (0:12) and Acos(awh)j from partons to measured 
jets. Since the correlation effects are stronger in jets that 
are closer together, the functions are parameterized in 
bins of cos(ai2) and cos(awb), respectively; one example 
for each function is shown in Fig. [2j 

The mw distribution after convolution with the func- 
tion Wa is shown in Fig. [3] The skewness is removed and 
the mean value agrees well with the parton-level distri- 
bution. 

The 20 integration variables (3 for each final-state par- 
ticle and the XBj for each initial state parton, assuming 
zero transverse momentum for the tt pair) in the expres- 
sion for the signal and background p.d.f.s (Eq. [2| are 
reduced to 16 by integrating over the 4-momcntum con- 
servation Dirac ^-function inherent in the expression for 
da s . The charged lepton 3- momentum integration and all 
but two of the jet angular integrations are made trivial 
by the Dirac ^-functions in the function W(k, x), leav- 
ing 7 integration variables. In P s , this is further reduced 
to 5 variables via a change of variables to the squared 
masses of the top quarks and by using the narrow-width 
approximation for the Breit-Wigner distributions of both 
top-quark decays in the tt matrix element. The integral 
is then evaluated using the VEGAS 24] adaptive Monte 
Carlo integration algorithm [6 , which uses importance 
sampling, which means that the sample points are con- 
centrated in the regions that make the largest contribu- 
tion to the integral. 

The treatment of Pj, is unchanged since the previous 
version of this analysis |BJ , except for the updated energy 
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FIG. 2: Examples of parameterization of the functions W\ 
and Wj h in the bins where 0.2 < cos (a 12 ) < 0.4 and 0.2 < 
cos (awb) < 0.4. The histograms show MC simulation events 
and the curves represent the parameterization. 



transfer function We- The integrand in the expression 
for Pf, is much more computationally intensive than for 
P s and a simplified Monte Carlo method of integration 
is employed, giving reasonable convergence with an ex- 
ecution time comparable to that of P s . The simplifica- 
tions used in this computation of Pf, include setting the 
function Wa to a Dirac <5-function for all angles, using 
a narrow width approximation for the W boson decay, 
and neglecting the Ajes dependence of the function We- 
Therefore, the value of P b for each event does not depend 
on the likelihood parameters m t and Ajes, while P s is a 
two-dimensional function of those parameters [6] . In this 
approximation, the product of the background p.d.f. nor- 
malization terms (corresponding to the variables &\, ■ Aj, 
in Eq.[2]) is set to a constant, whose value is chosen to op- 




FIG. 3: The mw distribution for measured angles from Fig.[T] 
is plotted (solid line) after convolution with the function Wa- 
For ease of comparison, the parton-level distribution (dotted 
line) is normalized so that the maxima of the two distributions 
are the same. 



timizc the statistical sensitivity of the method, effectively 
providing an appropriate relative normalization with re- 
spect to P s . 

The log-likelihood function is given as a sum over the 
578 events in the sample: 

ln£(fc; m t , A JES , i/ sig ) = (4) 

578 

y^ln [is sig P s (ki;m t , A JE s) + (1 - v sig )P b (ki)] . 

i=l 

It is calculated on a two-dimensional 31 x 17 grid in 
nit and Ajes, spanning 145 < m t < 205 GeV/c 2 and 
—4.8 < Ajes < 4.8, with a spacing between grid points 
of 2 GeV/c 2 in m t and 0.6 in Ajes- To optimize com- 
putational time, the bin size is chosen to be as large 
as possible without appreciably affecting the fit result. 
The third likelihood parameter, the signal fraction pa- 
rameter f s j g , is allowed to vary continuously (within the 
constraint < u S i S < 1), and the likelihood function is 
maximized with respect to v s - 1& at each point on the grid 
using the MINUIT program [25]. The resulting sur- 
face described on the grid is the profile log-likelihood, 
maximized for ^ s j g . The top-quark mass, m t , and the 
jet energy scale correction, Ajes, are measured by mak- 
ing a two-dimensional parabolic fit to the surface, con- 
sistent with the expectation for the likelihood function 
to be Gaussian near its maximum. The maximum of 
the parabola gives the measured m t and Ajes , while the 
measured u s i g is taken from its value at the grid point of 
maximum likelihood. The estimated one-a statistical un- 
certainty of the measurement is represented by the ellipse 



corresponding to a change in log-likelihood Aln£ = 0.5 
from the maximum of the fitted parabola. The values of 
m t and A JE s are anti-correlated (Fig. [4]). No correlation 
is observed between i/ s j g and mt or Ajes- 

The accuracy of the measured m t and Ajes, an d their 
uncertainties, are checked using ensembles of MC sim- 
ulated experiments, using the MC samples previously 
mentioned with the addition of 22 tt samples generated 
with values of m t between 161 and 185 GeV/c 2 . The 
numbers of events and those of the various backgrounds 
are Poisson fluctuated around the values shown in Ta- 
ble [Tl] Studies of the relationships between the known in- 
put simulation parameters and their corresponding mea- 
surements show no evidence of bias when a clean sample 
of MC simulated tt events is used, containing only lep- 
ton+jets events with correct jet-parton matching. How- 
ever, the presence of signal events with jets which are 
poorly or incorrectly matched to partons and events 
which do not match the decay hypothesis biases the like- 
lihood fit result and increases the pull width. The pres- 
ence of background events also biases the fit, due to the 
backgrounds that are not well described by P& and the 
approximations in P^. The bias is removed using a set 
of functions obtained from a fit to the MC simulation 
and parameterized in terms of the measured Ajes and 
v s i g [20]. This amounts to adding 1.1 GeV/c 2 to the mt 
value produced by the likelihood fit and multiplying the 
uncertainty by 1.26 so that the pull width is consistent 
with unity. The systematic uncertainty due to this mea- 



surement calibration is small, as shown in Table III 

Despite the reduction from the in situ Ajes calibra- 
tion, the remaining uncertainty from JES obtained by 
varying the parameters in JES [22] is among the largest 
systematic uncertainties of the measurement (Table III I. 



Other significant systematic uncertainties are mainly a 
result of assumptions made in the simulation of the events 
that are used in the tuning and calibration of the mea- 
surement method. In most cases, they are evaluated by 
varying different aspects of the MC simulation, such as 
signal MC generator (PYTHIA versus HERWIG [H]), 
color reconnection model tune (Apro versus ACRpro [27l - 
129] ). and parameters of initial and final state radiation 
(ISR and FSR) . A detailed description of the systematic 
effects has been published elsewhere [30.. The system- 
atic uncertainties for each effect are added in quadrature, 
resulting in a total estimated systematic uncertainty of 



1.3 GeV/c 2 (Table III) 



The measurement is made using the data sample of 
578 events, yielding 

m t = 172.4 ± 1.4 (stat+AjEs) ± 1.3 (sys) GeV/c 2 
mt = 172.4 ± 1.9 (total) GeV/c 2 , (5) 

with Ajes = 0.3 ± 0.3 (stat). The central value and 
the contour ellipses corresponding to the one-, two- and 
three-cr statistical confidence intervals of the measure- 
ment are illustrated in Fig. |4j The overall statistical 
uncertainty on the measured top-quark mass is labeled 



Systematic 


V / C ) 


MC Generator 


0.70 


Residual JES 


0.65 


Color Reconnection 


0.56 


fe-jet energy 


0.39 


Background 


0.45 


ISR and FSR 


0.23 


Multiple Hadron Interactions 


0.22 


PDFs 


0.13 


Lepton Energy 


0.12 


Measurement Calibration 


0.12 


Total 


1.31 



TABLE III: Contributions to the total expected systematic 
uncertainty. 
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FIG. 4: The measurement result and the contour ellipses of 
the parabolic fit corresponding to the one-, two- and three-cr 
confidence intervals for the statistical uncertainty on mt and 
Ajes. 



"stat+AjEs" because it includes the uncertainty on mt 
due to the statistical uncertainty on the measured Ajes j 
i.e. the uncertainty is given by half of the full width of 
the one-cr contour of Fig. [4] 

In conclusion, a precise measurement of the top-quark 
mass has been presented using CDF lepton+jets candi- 
date events corresponding to an integrated luminosity of 
3.2 fb . Using an improved matrix element method with 
an in situ jet energy calibration, the top quark mass is 
measured to be mt = 172.4 ± 1.9 GeV/ c 2 . 
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